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Abstract 

This paper considers control systems denned on Lie algebroids. After deriving basic controllability 
tests for general control systems, we specialize our discussion to the class of mechanical control systems on 
Lie algebroids. This class of systems includes mechanical systems subject to holonomic and nonholonomic 
constraints, mechanical systems with symmetry and mechanical systems evolving on semidirect products. 
We introduce the notions of linear connection, symmetric product and geodesically invariant subbundle 
on a Lie algebroid. We present appropriate tests for various notions of accessibility and controllability, 
and analyze the relation between the controllability properties of control systems related by a morphism 
of Lie algebroids. 

1 Introduction 

One of the basic problems in control theory is that of deciding the local controllability properties of a 
given system. Roughly speaking, the local controllability problem consists of finding appropriate conditions 
guaranteeing that the set of reachable states starting from an initial point is open, i.e., that the system can 
move locally in any direction. Deciding the controllability properties of a system is an a priori question that 
one needs to have addressed before being able to undertake other control problems such as motion planning 
and trajectory generation. The controllability problem has received a great deal of attention during the last 
decades (see El El El E| an d references therein). In particular, researchers have undertaken a thorough 
study of control systems with a rich geometric structure such as mechanical and homogeneous systems, and 
made use of their special properties to accomplish accurate modeling settings and sharp analysis results. 
Specific class of control problems include simple mechanical systems [2U], systems subject to nonholonomic 
constraints |3 03 El > systems invariant under the action of a Lie group of symmetries El El E] > systems 
enjoying special homogeneity properties [51 El EH , systems evolving on semidirect products (3U|, and more. 

One of the features which imposes a separate study for each class of systems is the lack of a unified 
framework. For instance, it is well known that a Lagrangian system invariant under the action of a Lie group 
of symmetries can be reduced to the quotient space induced by the action, but the reduced dynamics is not 
the one that corresponds to the reduced Lagrangian. Recent investigations have lead to a unifying geometric 
framework to overcome this drawback. It is precisely the underlying structure of Lie algebroid on the phase 
space what allows a unified treatment. This idea was introduced by Weinstein |H5| in order to define a 
Lagrangian formalism which is general enough to account for the various types of systems. A symplcctic 
formalism was later introduced for Lagrangian E] and Hamiltonian systems [21] ■ One of the advantages 
of the Lie algebroid formalism is the possibility of establishing appropriate maps (morphisms) between two 
systems that respect the structure of the phase space, and allow to relate their respective control properties. 
The underlying idea (in a similar way to the notion of abstraction 29 ) is that the property of interest will 
be easier to decide for one of the systems, and that by means of the morphism one will be able to infer the 
same knowledge for the other system. 
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In this paper we study the controllability problem for systems affine in the inputs evolving on a bundle E 
with an underlying Lie algebroid structure, r : E — ► M. We build on previous work on local controllability of 
general systems [22] and of mechanical systems 20_ to derive tests to check the accessibility and controllability 
properties of control systems evolving on a Lie algebroid. Throughout the paper, we pay special attention 
to what we term mechanical control systems on a Lie algebroid. This class of systems embraces a variety of 
different situations that can occur when analyzing mechanical systems, such as the ones mentioned above. 
Building on the notion of prolongation TE — > E of the Lie algebroid E — > M introduced in [25] , we develop 
all the necessary differential geometric tools enabling an intrinsic treatment of the second-order dynamics 
associated with mechanical systems. We focus on the set of reachable points in the base manifold M and 
in the bundle space E starting from states which belong to the zero section of E — > M. We make use of 
the geometric homogeneity properties of the controlled equations, which turn out to greatly simplify the 
accessibility and controllability computations. We carefully describe the relation between the controllability 
properties of control systems that are related by a morphism of Lie algebroids. As a result of the generality 
of the approach, we are able to present in a unified way previous work on the configuration accessibility and 
controllability properties of simple mechanical control systems E0 EOI ES] (see also 0] for a comprehensive 
overview). Regarding systems evolving on semidirect products, the application of the Lie algebroid approach 
renders novel tests which are valid in slightly more general settings than the ones considered in [20] ■ We 
also extend notions such as fiber controllability to what we call controllability with regards to a manifold and 
develop conditions to check this property. 

In the course of the preparation of this manuscript, we came across the recent research effort |34|. This ref- 
erence, which is close in spirit to this work, analyzes the controllability properties of so-called "1-homogeneous 
control systems" evolving on a vector bundle. However, it deals with vector fields with values in the tangent 
bundle of the vector bundle, as opposed to deal with the formalism of Lie algebroids and their prolongations. 
This choice of a higher-dimensional phase space makes necessary to resort to additional geometric tools such 
as Ehresmann connections in order to describe the structure of the accessibility algebra. We think that the 
Lie algebroid approach accommodates the same level of generality, while enabling in general a more concise 
treatment of the controlled dynamics. 

The paper is organized as follows. In Section [2] we present some basic facts on Lie algebroids. We 
also discuss in detail the notion of linear connection on a Lie algebroid, including the generalization of 
the Levi-Civita connection and the constrained connection. In Section [2] we introduce the prolongation 
of a Lie algebroid and develop the differential geometry of horizontal sections, homogeneity, SODE sections 
and geodesically invariant subbundles. In Section 0] we study nonlinear affine control systems whose drift 
and input vector fields are associated with some sections of a Lie algebroid. This apparent restriction is 
not such, since most physical systems can be casted into this form. We formulate the conditions for local 
accessibility and controllability in terms of Lie brackets of sections, and we study the effect of a morphism of 
Lie algebroids on these properties. In Section we introduce the class of mechanical control systems defined 
on a Lie algebroid. We show that the notion of affine connection control system can be generalized to the 
setting of Lie algebroids, thus providing a general framework to study the controllability properties of these 
systems. We introduce the notions of local base controllability and controllability with regard to a manifold, 
and we obtain computable sufficient conditions to check them. We also study the effect of morphisms of Lie 
algebroids in simplifying the controllability analysis. These results are later applied in Section [S] to simple 
mechanical systems defined on a manifold, simple mechanical systems with symmetry and systems defined 
on semidirect products and orbits of group actions. Section presents some concluding remarks. We have 
gathered in an appendix some basic notions for control systems defined on manifolds. A final remark is that 
the summation convention over repeated indexes is understood throughout the paper. 

2 Preliminaries on Lie algebroids 

In this section we introduce some known notions and develop new concepts concerning Lie algebroids that are 
necessary for the further developments. We refer the reader to [23 [52] for thorough studies of Lie groupoids, 
Lie algebroids and their role in differential geometry. Let M be an n-dimensional manifold and let r : E — ► M 
be a vector bundle with ^-dimensional fibers. A structure of Lie algebroid on E is given by a Lie algebra 
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structure on the C 00 (M)-module of sections of the bundle, (Sec(E), [•,•]), together with a homomorphism 
p: E — > TM of vector bundles (called the anchor map) satisfying the compatibility condition 

[ai,Fa 2 ] = F[a 1 ,a 2 ] + {p{ax)F)a^ . 

Here F is a smooth function on AI, <ti, a 2 are sections of E and we have denoted by p(cr) the vector field 
on M given by p(o~)(m) — p(o~(m)). The homomorphism p is called the anchor map. From the compatibility 
condition and the Jacobi identity, it follows that the map a i— ► p(a) is a Lie algebra homomorphism from 
Sec{E) to X(M). 

It is convenient to think of a Lie algebroid r: E — > M as a substitute of the tangent bundle of M. In 
this way, one regards an element a of E as a generalized velocity, and the actual velocity v is obtained when 
applying the anchor to a, i.e., v = p(a). 

The image of the anchor map, p(E), defines an integrable smooth generalized distribution on M. There- 
fore, M is foliated by the integral leaves of p(E), which are called the leaves of the Lie algebroid. A curve 
a : [ta,ti] — > E is said to be admissible if m(t) = p(a(t)), where m(t) = r(a(t)), t £ [ta,ti\. It follows that 
a(t) is admissible if and only if the curve m(t) lies on a leaf of the Lie algebroid, and that two points are in 
the same leaf if and only if they are connected by (the base curve of) an admissible curve. 

A Lie algebroid is said to be transitive if it has only one leaf, which is obviously equal to M. It is easy 
to see that E is transitive if and only if p is surjective. If E is not transitive, then the restriction of the 
Lie algebroid to a leaf L C M, E\l — > L is transitive. In the latter case, one can show that E\l, and hence 
ker p, has constant dimension. We will say that a Lie algebroid is locally transitive at a point to 6 M if 
p m : E m — > T rn M is surjective. In this way, to is contained in a leaf of maximal dimension. 

Given a local basis {e Q }^ =1 of sections of E defined on an open set V C M, we can write a = y Q e Q (r(a)) 
for any a £ E such that r(a) G V. If (x 4 ), i = 1, . . . , n are local coordinates in the base M defined on V, 
we have local coordinates (a;*, y a ), i = 1, . . . , n, a = 1, . . . , £ in E. The anchor map and the Lie bracket are 
then determined by the local functions p l a and Cg 7 on M (called the structure functions of the Lie algebroid) 
defined by 
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The structure functions satisfy the following relations 

cyclic(a,/3,7) 



l PI , ftfi fiv 



(2.1) 



where the summation over repeated indexes is understood. Equations l|2. 1 1> are usually called the structure 
equations of the Lie algebroid. Finally, the Lie bracket of two sections of E can be expressed in terms of the 
basis {e a Ya=i as 

k »] - (°"P*% ~ ^7 15 + ^> V ) ^ ■ (2 ' 2) 

If y is a family of sections of E, we will denote by Lie(y) the distribution obtained by closing (the distribution 
defined by) y under the Lie bracket. 



2.1 Admissible maps and morphisms of Lie algebroids 

Let r: E — > M and r: E — * M be two Lie algebroids with associated anchor maps p: J5 — > TM and 
p: i? — > TM. A bundle map \I/ : £7 — > £7 is said to be admissible if op = p o vp. Equivalently W is 
admissible if and only if it maps admissible curves into admissible curves. Indeed, if a(t) is admissible on E 
and projects to m(t), then a(t) = *£(a(t)) projects to m(t) — ip(m(t)) and it is admissible, since 

p(a(t)) = p(*(a(t)) = T^j(p(a(t)) = T^(m(t)) = m(t). 
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Denoting by ip: M — > M the map on the base, one has the following commutative diagram 




A map \P : E — > E is a morphism of Lie algebroids if it is admissible and preserves the Lie algebra structure 
of the algebroids that is, for any a and rj sections of E such that there exist some sections {0}f = i of E 
and some functions Fi , Gi , ^ = 1, . . . , p on M with 



(=1 

then, the image of the Lie bracket of a and 77 under 4" is 

p 



z=i 



2=1 



Given local basis {e a } e a=1 and {e a } e a=1 of sections of 15 and E, respectively, a bundle map ^ can be written 
$(e„) = ^^ep f° r certain local functions on M. Then, one can check that W is a morphism if and only 
if 



' 7 a5 



dx l 



(2.3) 



Notice that if a, r\ are VP-related to sections a, r\ e Sec(-E), i.e., $ o a — a o -0 and $ o t/ = ij o ^i, then the 
Lie bracket [<r, 77] is ^-related to the Lie bracket [<r, 77], \& o [cr, 77] = [er, 77] o 0. 



2.2 Linear connections 

Here we briefly present the notion of ^-connection on a vector bundle (cf. [TT], see also 0^]), and discuss 
some related objects. 

Definition 2.1. Let t : E —* M be a Lie algebroid. A linear E-connection on a vector bundle ir: P — > M is 
a M.-bilinear map V: Sec(E) x Sec(P) — > Sec(P) such that 

V F<J a = FV a a and V a {Fa) = {p(a)F)a + FV a a 

for any function F G C°°{M), section a of E and section a of P. 

Throughout the paper, we will restrict our attention to the case P — E, and by a connection on E wc 
will understand a linear ^-connection on t : E — > M. Given a local basis {e a } e a=1 of sections of E, the local 
expression of the covariant derivative is 

i drj 
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The terms r^g are called the connection coefficients. As in the study of tangent bundle geometry, the 
skew-symmetric part of the connection defines the so-called torsion tensor, 

T(a, rj) = V CT 77 - V„ct - [a, n) , 

and the symmetric part of the connection determines what we call the symmetric product, 

(a :rj) = S7 a n + V, ; cr . 

The local expression of the symmetric product is (compare with the expression for the Lie bracket (|2.2[l l 

/ f) n a f> n ct \ 

{a ■ v) = \ alp "-i? + 7fp "d^ + ^W) e « ' where s ^ = r ^ + r k • ( 2 - 4 ) 

In particular, notice that (e Q : ep) = S'2/3 e 7- Similarly as with the involutive closure, if y is a family of 
sections of E, we will denote by Sym(^) the distribution obtained by closing (the distribution defined by) V 
under the symmetric product. 

Since the covariant derivative is C°°(Af)-linear in the first argument, it is possible to define the derivative 
of a section a G Sec(E) with respect to an element a G E m by simply putting 

V a fT = V,,Ct(to) , 

where w G Sec(-E') is any section such that rj(m) = a. Moreover, the covariant derivative allows us to take the 
derivative of sections along maps and, as a particular case, of sections along curves. If we have a morphism 
of Lie algebroids <I>: F — > E over the map <p: N — > M, then we can define the derivative of a section of E 
along ip as follows. 

Definition 2.2. Let a : N — > E be a section of E along p>, i.e., o~(n) G £w„), n G N . Then a can 
be written in the form a = X)f=i^ 7 i(0 <p)> f or some sections {£i,...,Cp} C Sec(£7) and some functions 
{Fx, . . . , F p } C C°°(N). The derivative of a along ip is defined by 

p 

V6^ = £[0 , ( 6 ) f »)Ci(v(n))+Ji(n)V#(6)Cl] . beF. 
1=1 

Remark 2.3. Within this framework, one can consider time-dependent sections of E as follows: take the 
morphism $: TR x E — > E, <&(t,a) = a over the map tpiM.xM — > M, p>(t,m) — m. The Lie algebroid 
structure on TR x E is the direct product structure, that is, the anchor is Ptrxe(t-^, a) = rj^ + p(a), and 
the bracket of projectable sections (on both factors) is the sum of the brackets on TR and E. 

When studying mechanical control systems related by a morphism of Lie algebroids, we will resort to the 
following notion concerning the interplay between maps and linear connections. 

Definition 2.4. Let V and V be connections on E and E, respectively, and let ^ be a bundle map from E 
to E. We say that "J maps the connection V to the connection V if 

*° (V*V) = V ct (*ot7). 

In coordinates this condition is equivalent to 

ami 3 a 

^r 7 x = p 1 — s - + Fn ^ e m . (2.5) 

Geodesies 

Consider the following situation: let a : 1 1— > E be an admissible curve, and let b : 1 1— > E be a curve in E, both 
of them projecting by r onto the same base curve in M, r(a(t)) — m(t) — r(6(i)). Take the Lie algebroid 
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structure TR — > R and consider the morphism $ : TM -> £, = ta(t) over 99 : R — » M, p(t) = m(t). 

Then one can define the derivative of b(t) along a(t) as V_d_6(i). In the literature, this derivative is usually 
denoted by V a (t)6(t). In local coordinates, this reads 



V„(t)&(t) 



e 7 (m(i)). 



Definition 2.5. Let t : E —> M be a Lie algebroid and V a connection on E. An admissible curve a : 1 1— > 
is said fo 6e a geodesic of V ifV a ^a(t) = 0. 



In local coordinates, the conditions for being a geodesic reads 



dt 



(2.6) 



The Levi-Civita connection 



Let S : E Km E — > R be a bundle metric on a Lie algebroid t : E ^ M. In a parallel way to the situation 
in tangent bundle geometry, one can see that there is a canonical connection on E associated with 9- The 
proof is analogous and will be omitted. 

Proposition 2.6. Given a bundle metric 9 on E , there is a unique connection on E which is torsion-less 
and metric with respect to 9- The connection is determined by the formula 

29(V CT r?, C) - p(o-)9(v, + PivUa, Q - P«)S(v, a) + 9(a, [£ rj\) + Sfa, [C, <A) ~ S(C, [»», *]) , 
/or 17,77, C G Sec(iJ). 

Denoting by {e Q }^ =1 a local basis of sections of E, and by {e a } e a=1 its dual basis, the bundle metric can 
be locally written as 9 = 9 a /3 e a ® e@ . The connection coefficients of V s are 

where (9^) is the inverse matrix of (Sap), and [a, f3; 7] is a shorthand notation for 

Associated with the bundle metric 9> one has the musical isomorphisms 

b s :E^E*, (b s (a), b) = 9(a, 6), fl g : E* ^ E , j) g (0) = ^(0) , 

where £J* — > M denotes the dual bundle of E. Given a function V on M, the gradient of V, gradg V, is the 
section of E defined by grad g V = $${p*dV). 



The constrained connection 



Here, we introduce the notion of constrained connection on a Lie algebroid, which generalizes the concept 
of nonholonomic connection ^JGSl- This notion will be later useful to model control systems subject to 
nonholonomic constraints, following the developments in ^ I19| . Given an arbitrary connection and an 
arbitrary subbundle of a Lie algebroid, one can define a new connection which enjoys special properties with 
respect to the subbundle. This connection is determined by the choice of a projector. Let D be a subbundle 
of E and let P be a projector onto D, P : E — > D. We denote by Q the complementary projector of P, 
Q = I — P, and by D c the complementary subbundle D c = Im(Q). In this way, one has D © D c = E. 

Definition 2.7. Given a connection V on E and a projector map P : E — » D, the constrained connection 
is the connection V on E defined by 

V CT ry = P{V a n) + V CT (Qr?) , tr,ij6 Sec(£) . 
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Some interesting properties of the constrained connection V are the following. Their proof is straightfor- 
ward and will be omitted for brevity. 

Proposition 2.8. The following properties of the constrained connection V hold: 

(i) The connection V restricts to D, i.e., V a r\ £ D for any 77 S Sec(-D) and a G Sec(iS). 
(ii) The symmetric product (■ : ■}" associated with V is given by 

• (a : 77}"= P((a : 77}) for a,r) G Sec(L>). 

• (o- : 77}"= P((<7 : 77}) + (cr : t;} /or 0", 77 G Sec(L> c ). 

• (o- : 77}-= P((o- : 77}) + V ff 77 for a G Sec(D) and 77 G Sec(D c ). 

(Hi) The torsion tensor T o/V is given by 

. f(a, r,) = P(T(a, 77)) - Q([a, 77]) /or a, 77 G Sec(D). 

• f (a, 77) = P(T(<7, 7/)) + V CT 7/ - V n< 7 - Qfljr, 77]) /or a, 77 G Sec(I? c ). 

• f (a, 77) = P(T(<7, 77)) + - Q([ct, 77]) /or a G Sec(L>) and 77 G Sec(D c ). 

Proposition 2.9. Let r : E — > M, t : E —> M be two Lie algebroids, with projectors P : E D and 
P : E — > £>. Lei V and V be connections on E and E, respectively. Assume that a morphism of Lie 
algebroids * : E -^E maps V onto V. J/*oP = Pof (equivalently C Z> and *(L» C ) C £> C J, i/ien * 

maps V onto V. 

Proo/. Since Q = I - P, it follows that f oQ = Qof. Therefore, for all 77 G Sec(E) and b e E 
*(V 6 77) = *(P(V 6 77)) + *(V 6 (Q ? /)) = P(*(V 6 77)) + V 6 (* o (Qn)) 
= P(V b (* o 77)) + V b (Q(* o 77)) = V h (* o ?/ ) . 

□ 

3 The prolongation of a Lie algebroid 

Here we briefly review the notion of the prolongation of a Lie algebroid. For further details, see |25|. Given 
a Lie algebroid E, the underlying motivation behind the introduction of the prolongation of E is that of 
formulating second-order dynamical systems on E. Thinking of E as a substitute of the tangent bundle of 
M, the tangent bundle of E is not the appropriate space to describe second-order dynamics on E. This is 
clear if we note that the projection to M of a vector tangent to E is a vector tangent to M, and what one 
would like instead is an element of E, the new tangent bundle of M. 

A space which takes into account this restriction is the i?-tangent bundle of E, also called the prolongation 
of E, which we denote by TE. This Lie algebroid is defined as the vector bundle n : TE — > E whose fiber 
at a point a G E m is the vector space 

T a E = { (6, v) G E m x T a E \ p(b) = T a r{v) } . 

Note that if the fibers of E are ^-dimensional, then the fibers of TE are 2^-dimensional. We will use the 
redundant notation (a, 6, v) to denote the element (6, v) G T a E. 

The anchor of TE is the map p 1 : TE — > TE, defined by p 1 (a, b, v) = v. We also consider the map 
Tt: TE — > E defined by Tr(a, b, v) = b. The Lie bracket associated with TE is defined as follows in terms of 
projectable sections. A section Z of TE is projectable if there exists a section a of E such that Tt o Z = a or. 
Equivalently, a section Z is projectable if and only if it is of the form Z(a) = (a,a(T(a)),X(a)), for some 
section a of E and some vector field X on E. The Lie bracket of two projectable sections Zi and Z 2 is then 
given by 

[Z u Z 2 ](a) = (a, [a 1 ,a 2 ]{m), [X 1 ,X 2 ](a)) , aeE. 
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It is easy to see that [Z\, Z 2 ](a) is an element of TE for every a G E. Since any section of TE can be locally 
written as a linear combination of projectable sections, the definition of the Lie bracket for sections of TE 
follows. 

Given local coordinates (x l ,y a ) associated with a basis {e Q } of sections of E, we can define a local basis 
{X a , V a Ya=i °f sections of TE by 



X a (a) = {a,e a (T{a)),p l a -^- J and V a (a) = (a, 0, ^ J , a = l,...,£. (3.1) 

If (a, 6, f ) is an element of TE, with 6 = z a e a and w = p l a z a J~ + v a -^^, then we can write 

(a,6,w) = z Q <Y Q (a)+w Q V Q (a). 

The Lie brackets of the elements of the basis are 

[X ol ,X f3 ] = Cl f) X 1 , [X a ,Vp}=0, [V a ,Vp}=0. 

Finally notice that the anchor map p 1 applied to a section Z = Z a X a + V a V a of TE defines a vector field 
on E whose coordinate expression is 

The following notion shows that morphisms of Lie algebroids can also be prolonged. 

Definition 3.1. Let ^ : E — > E be a an admissible map between two Lie algebroids r : E — > M and t : E — > 
M. The prolongation of * is the mapping T*: TE — » TE defined by T^(a,b,v) = (*(a), T a V(v)). 

It is not difficult to see that if is a morphism of Lie algebroids, then its prolongation is also a morphism of 
Lie algebroids. Given local basis {e Q }^ =1 and {e a }^ =1 of sections of E and E, respectively, if $(e tt ) = ^fe^, 
then the action of the prolongation of W on TE is determined by 

T^(X a (a)) = *g(m)^(*(o)) + p l a (m)^{m)a^V^{a)) , a G E m , 
r*(Va(a)) = *g(m)V /J (*(o)), «e£ m . 
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3.1 Vertical and horizontal sections 

An element (a,b,v) G T2£ is said to be vertical if it is of the form (a,0,v), with w a vertical vector tangent 
to E at a. It follows that the vertical space of TE at the point a G E m , which we denote by Ver a (TE), can 
be identified with E m by (a slight modification of) the usual vertical lifting map: 

b G E m 1 — ► (a, 0, &£) G Vcr a (T£) , 

where b v a G T a E is the tangent vector to the curve t ^ a + tb at t = 0. If cr is a section of -E 1 , then the section 
g v of TE defined by <J V (a) = (a,0,a(m)^) will be called the vertical lift of a. Vertical elements are linear 
combinations of {V Q }^ =1 . Specifically, if the section a of E has the local expression a = a a e ai then o v is of 
the form 

a v = a a V a . 

Consider the zero-section of r : E — > M, that is 

0m:M^E, M (m)=0 m . 

This section is a canonical embedding of M into E. Consequently, we can regard TM as a subspace of TE. 
Now, define the horizontal space along 0m, 

Hor m (T£) = { (0 m , b,v)e% m E \ ve T m M CT 0m E}, meM. 

Note that Ti(Hov(TE)) = Im(OM), i-e. Hor(TE) is only defined along points of E that belong to the 
zero-section M . Moreover, Hor m (TE) is a vector subspace of T 0m E. The following result shows that the 
horizontal space is complementary to the vertical space along the zero-section 0m- 
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Lemma 3.2. Along the zero-section of r : E — > M, we have the following direct sum decomposition 

% M E = Hor(TE) © Vct 0m (TE) . 

Proof. Consider the map "H: E — > % AI E given by J£(&) = (0 T (6), 6, T0m(p(6)))- The image of IK is precisely 
the horizontal space, Hor m (2"E) = 3i(E m ), m G M. Note that 3~C is a splitting of the short exact sequence 

»- £ >■ %u E " E *~ , 

where V(b) = (0 r ( b ), 0, 0$ ). As a consequence, we have T Qm E = 0i(E m ) © Ver 0m (TE), m G M. The 
restriction of 7r to Hor(T_E) is an isomorphism 7r : Hor(T_E) — * E 1 , whose inverse map is Jf . For (0 m , 6, u) G 
7o m i?, the decomposition is given by (0 m ,b,v) = (0 m , 6, T(Qjvfor)(u))+(0 m , m , v-T(0 m °t)(v)) G Hor(T_E) + 
Ver 0M (TS). □ 

3.2 Homogeneity 

A property that plays an important role in our later analysis is that of homogeneity. Consider the section 
A of TE defined by A (a) = (a, 0, a„), a G The section A is called the Liouville section of TE. In 
coordinates, we have 

A = y a V a and p 1 (A) = y«^. 
A function _F G C°° (E) is said to be homogeneous of degree s G Z if 

&pi(A) F = sF ■ 

where £ stands for the Lie derivative operator 1 . In a local chart, a homogeneous function of degree s > 
is a homogeneous polynomial in {y a }a=i of degree s with arbitrary functions of (x 1 )™^ as coefficients. 
Consequently, homogeneous functions of degree are (pullbacks of) functions on the base M and there 
are not (smooth) non-trivial functions homogeneous of degree s < — 1. A section Z of TE is said to be 
homogeneous of degree s G Z if 

[A, Z] = sZ. 

We denote by CP S the set of homogeneous sections of TE of degree s. The following result describes the basic 
properties concerning homogeneous sections. The proof is omitted for brevity. 

Lemma 3.3. Let r, s G Z and let Z be a section ofTE. Then 
(i) [9 s ,9 r ] C y s+r , and 7 S = {0} ifs<2, 

(ii) Z G CP i if and only if there exists a section a of E such that Z = o~ v , 

(Hi) Z G To if an d only if Z is a projectable section, 

(iv) X a G To and V a G for a = 1, . . .,£. Moreover, if Z = Z a X a + V a V a is the local expression of 

Z , then Z G CP S if and only if the functions Z a are homogeneous of degree s and the functions V a are 
homogeneous of degree s + 1. 

Note that for all Z G CP S , s > 1, we have Z(0 m ) = Oo m , m G M, that is, the homogeneous sections of TE 
of degree greater or equal than 1 vanish at the zero-section of E. 

1 Alternatively, one can see that homogeneous functions of degree s verify d&F = sF. (See 1211 for the precise definition of 
the derivative operator d&). 
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3.3 SODE sections 



The dynamics of a mechanical system evolving on a certain configuration manifold is described by means of 
a vector field on the tangent bundle of the manifold which is a second order differential equation. Likewise, 
we will need the notion of SODE section of the prolongation to describe the behavior of mechanical systems 
evolving on Lie algebroids. In this section, we introduce this concept and discuss several geometric properties, 
resembling those of second order differential equations on manifolds. 

Definition 3.4. A section T ofTE is a second-order differential equation ( SODeJ section on the Lie algebroid 
E ifTroT = Id E - 

A vector v 6 T a E is called admissible if Tp(v) — p(a). Note that a curve in E is admissible if and only if 
its tangent vectors are admissible. We denote by Adm(iS) the set of all admissible tangent vectors. Notice 
that v is admissible if and only if (a, a, v) £ TE. Therefore we can identify Adm(_E) with the subset of TE 
formed by all elements of that form, 

Adm(£) = {z e TE | n(z) = Tr(z)} . 

Equivalently, a SODE section can be defined as a section of TE that takes values in Adm(E), i.e., T(a) — 
(a, a, X(a)), a 6 E, where X is a vector field on E verifying p(a) — T a r(X(a)). If L is a SODE section, then 
it has the local expression L = y a X a + F a (x, y)V a and its associated vector field is of the form 

The integral curves of this vector field satisfy the differential equations 

x^piv", y a = F a (x,y). 

In particular, we are specially interested in homogeneous SODE sections with degree 1. We will refer to such 
SODE sections as sprays. Locally, a spray is such that the functions F a are homogeneous with degree 2, 
F a (x, y) = —^Sp 1 {x)y^ i y 1 for some symmetric coefficients Sp . The sprays are in one to one correspondence 
with torsion- less £J-covariant derivatives (cf. Section |2~2"|) . as we show in the following. 

Let V be a connection on E. The geodesic equations (|2.(j|) . together with the admissibility condition, 
correspond to the differential equations of the integral curves of a SODE section Ly of TE, which is locally 
given by 

r v = ^A" a -i(r| 7 + r^) yVv a . 

From the coordinate expression, we easily see that Fy is homogeneous with degree one, i.e., it is a spray. 
Moreover, notice that Ty is determined by the symmetric product associated with V (cf. eq. 12.4| l). 

We now show how a spray T and a (2,l)-tensor field T determine a connection on E in a unique way. As 
an intermediate step, we first define the symmetric product associated with a spray. Given a, n € Sec(E), 
consider the section of TE, [a v , [r, r) v ]]. From the properties of the Lie bracket, we deduce that this section 
is (smooth and) homogeneous with degree —1. From Lemma ISTST ii'). it follows that [a v , [T, n v ]} is the vertical 
lift of a section of E, which we denote by (a : rj)r- Therefore, 

{a:n) v = [o- v ,[T,n v ]]. 

From the Jacobi identity, one can deduce that the operation (a : rj) is symmetric. Moreover, if / is a function 
on M, we have (a : frj) = p(a)fr) + f(cr : 77), since 

(a : fnY = [K,r], (/„)"] = p l {[a\T))frf + f[a v , [T, n v }) = (p(a)fr, + f(a : n)) v . 

Definition 3.5. Given a spray T on E, (■ : -)r is the symmetric product associated with T. 
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Taking a local basis of sections of E, one can compute 

[r, a v ] = [y a X a - \s%^yiV a yVp] = -a a X a + y a + ST^tA V 7 , (3.3) 

[V v , [r, a v ]] = [e'fijk + + ^Vj V a . (3.4) 

From the local expression of the symmetric product, one can see that the symmetric product determines and 
is determined by the spray T. 

Proposition 3.6. Given a spray T and a skew- symmetric (2,1) tensor field T on E, there exists a unique 
connection V r ' T : Sec(-E) x Sec(-E) — > Sec(-E) on E such that its associated spray is V and its torsion is T. 
This connection is given by 

V T a T V = \ {[a, rj\ + T{a, rj)) + \{a:r 1 ) T . (3.5) 

Proof. From the properties of the symmetric product and the Lie bracket, it follows that V r,T defined in l|3.5[) 
is a connection on E. Its torsion is given by 

V T a T V - W r n T a - [a, r,} = \ ([a, rj\ + T(a, r,)) -\{%a} + T(r h a)) - [a, v ] = T(a, rf). 

By definition, the spray associated with V r ' T is the SODE section Z whose projection p l {Z) is determined 
by the differential equation V a ^a(f) = 0, for admissible curves a : K — > E, t i— > a(t). A simple coordinate 

calculation shows that these equations can be locally written as y a + ^S^y^y 1 — 0, which concludes the 
proof. □ 

The connection coefficients of V r ' T are given by 

L f3~f- j V ^ Pi P~< 

Note also that the symmetric product associated with T precisely corresponds to the symmetric product 
defined by Vr t- 

(a:r l ) T = V v a T V + V^ T a. 

Proposition 3.7. Let r : E — > M and t : E — > M be two Lie algebroids. Let V and V be connections on E 
and E, respectively. A morphism ^ maps V onto V if and only if ^ maps the associated spray T into the 
associated spray T and maps the torsion tensor T into the torsion tensor T. 

Proof. We prove in coordinates the above statement for torsion-less connections, from where the general 
result follows easily. If we take a SODE T — y a X a + F a V a , then 



r*(r(o)) = a a x a 

This last expression is equal to T(a) = a a X a + F~(a)V a if and only if 



Xy , a G Er n 



In the case of two sprays, F 1 — —^S2 l gy a y^ and F 7 = —^S 1 a py a y 13 , so that the above equation reads 



Since \1/ is a morphism, it verifies equation l|2.3|l . Summing both expressions, we finally obtain equation (|2.5|) . 
as claimed. □ 
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3.4 Geodesically invariant subbundles 



Here we introduce the notion of geodesically invariant subbundles (which is a generalization of the concept 
of geodesically invariant distributions |18|L and establish its relation with the symmetric product associated 
with the connection. 

Definition 3.8. Let V be a connection on E. A subbundle D C E is geodesically invariant for V if every 
geodesic a: [io,ii] — > E such that a(to) G D verifies that a(t) G D for all t G [to,ii]. 

Proposition 3.9. Let V be a connection on E. A subbundle D C E is geodesically invariant for V if and 
only if it is invariant under the symmetric product, i.e. (a : rf) G Sec(D) for every cr, 7/ G Sec(Z?). 

Proof. By definition, the subbundle D is geodesically invariant if and only if the spray of the connection 

IV is tangent to D, i.e. p 1 (T)\d G TD. The latter is equivalent to the condition (£j p i (r)0) \ D — for 
every function (j> on E such that 4>m = 0. Since D is linear, the constraint functions <p are of the- form 
4> = fi, with p G Sec(-D°) and p(a) — (p m ,a), a G E m . Therefore, D is geodesically invariant if and only 
if (X p i( r )/i) o a — for all p G Sec(Z?°) and a G Sec(-D). Let us see that (L p i^p) ° °~ = —^(fa^o-a). In 
coordinates 

Taking the restriction to the image of a, and writing explicitly the symmetric parts, we have 

V(r)M o a = - [Pp-Qj + P a -gj ~ PiSl ja a". 

On the other hand 
and thus 

where we have used the fact that + = 52^. Therefore (£ p i(r)A) ° a — 5(V ff jti, it), and the result 
follows by taking into account that (V CT jU, cr) = p(a)(/j,, a) — (/i, Vo-ct) = — (/i, V CT cr), because (/i, cr) = 0. □ 

Given a connection V on E and a projector map P : E D, consider the constrained connection 

V introduced in Section 12.21 Using the above result and Proposition I2.8f ii) , one can deduce that D is 
geodesically invariant for V. 



4 General control systems on Lie algebroids 

In this section we present the notion of a control system on a Lie algebroid. We introduce the concept of 
accessibility subbundle in the Lie algebroid and provide basic tests for controllability, building on the known 
results for control systems defined on manifolds. Finally, we study the controllability properties of control 
systems related by means of a morphism of Lie algebroids. 

Consider a Lie algebroid r : E — > M, with anchor map p : E — > TM. Let cr, r/i, . . . , rjk be sections of E. 
A control problem on the Lie algebroid E — > M with drift section a and input sections 771, . . . ,T)k is defined 
by the following equation on M, 

k 

m(t) = p(a{m{t))) + £ Ui {t) m {m{t))) , (4.1) 
i=l 

where u = (u\,...,Uk) G U, and U is an open set of R fc containing 0. The function t u(t) = 
(ui(t), . . . ,u m (t)) belongs to a certain class of functions of time, denoted by U, called the set of admissible 
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controls. For our purposes, we may restrict the admissible controls to be the piecewise constant functions 
with values in U. Notice that the trajectories of the control system are admissible curves of the Lie algebroid, 
and therefore they must lie on a leaf of E. It then follows that if E is not transitive, then there are points 
that cannot be connected by solutions of any control system defined on such a Lie algebroid. In particular, 
the system 1(4.1)) cannot be locally accessible at points to G M where p is not surjective. Since the emphasis 
here is put on the controllability analysis, without loss of generality we will restrict our attention to locally 
transitive Lie algebroids. 

Denoting by / = p(a) and gi = p(r)i), we can rewrite the system (|4.1(l as 

k 

m(t) = f(m(t)) + ui(t) gi {m(t)) , (4.2) 
»=i 

which is a standard nonlinear control system on M affine in the inputs. Here we make use of the additional 
geometric structure provided by the Lie algebroid in order to carry over the analysis of the controllability 
properties of the control system 1)4.1(1 . We refer to [2E] for a comprehensive discussion of the notions of 
reachable sets, accessibility algebra and computable accessibility tests. A short list of definitions is provided 
in the appendix for reference. 

Definition 4.1. The accessibility algebra D of the control system 1)4. If) in the Lie algebroid is the smallest 
subalgebra of Sec(E) that contains the sections a, 771, . . . , 

Using the Jacobi identity, one can deduce that any element of accessibility algebra D is a linear combi- 
nation of repeated Lie brackets of sections of the form 

[0, [...,&,&] •••]]], 

where d G {a, 771, . . . , 7^}, 1 < i < I and I G N. 

Definition 4.2. The accessibility subbundle in the Lie algebroid, denoted by Lie({cr, 771, . . . , Tjk}), is the 
vector subbundle of E generated by the accessibility algebra D, 

Lie({cr, 771, . . . , rjk}) = span {C(to) | C section of E in D} , to G M . 

If the dimension of Lie({<7, 771, . . . ,7/fe}) is constant, then Lie({er, 771, . . . ,77k}) is the smallest Lie subalge- 
broid of E that has {a, 7/1, . . . ,r/k} as sections. In the following result, we establish accessibility tests for 
control systems of the form ()4.1|) . 

Theorem 4.3. Consider the system 114.1(1 . Let to G M and assume the Lie algebroid E is locally transitive 
at 771. Then, Lie({er, 771, ... , rjk}){m) + ker p(m) = E rn implies that the system is locally accessible from m. 

Proof. We first prove that every element X of the accessibility algebra 6 of 14. 2|) can be written as X = p(C), 
with ( G B (see the appendix for the precise definition of the accessibility algebra C). Take an element 
of 6 of the form [X h [. . . , [X 2 ,X 1 ] . ..]]], with X t G {/ = p(o-),gi = p(m), ■■■,9k = p(Vk)}- Denote 

Xi = p(Ci)i with d G {<r, 771, . . . , 77fe}. Since p is a Lie algebra homomorphism, we have that 

[X U [AVi, [.-. , [X 2 ,X 1 ] ...}}} =p([Q, [. . . , [C 2) &] ■ ••]]]) • 

Since p is linear, we conclude the accessibility subbundle in the Lie algebroid is mapped by p onto the acces- 
sibility distribution C . Now, we show C(m) = T rn M if and only if Lie({<7, 771, ... , rji~})(m) + ker p(m) — E m . 
Assume first that C(m) = T m M. Let e m G E m . Consider p(e) G T m M — C(m) — p(Lie({er, 771, ... , r]k})(m)). 
There exists a G Lie({cr, 771, ... , 7^}) (to) such that p(e) — p(a). Then, e — a G ker p(m) and e = a + (e — a). 
Therefore E m = Lie({cr, 771, . . . , rjk}) (m) + ker p(m). In addition, p(E m ) = p(Lie({er, 771, ... , ?7a;})(to)) = T m M. 
The other implication is trivial. Finally, the result follows from Chow's theorem [2H1- D 

Remark 4.4. If the dimension of Lie({<7, 771, . . . , 77^}) is constant, then the above theorem expresses the 
following fact: if the Lie subalgebroid Lie({er, 771, . . . , 77^}) is locally transitive at m, then the system (|4.1(l is 
locally accessible from to. 
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In practice, the most interesting property to establish is controllability. In what follows, we provide a 
controllability test for systems of the form 1)4.111 by adapting the notions of good and bad Lie brackets of 
vector fields to the setting of Lie algebroids. Let {Xo, Xi, . . . , X/.} be a set of sections of the Lie algebroid 
E. The degree of an iterated Lie bracket B of elements in {Xq, Xi, . . . , Xk} is the number of occurrences of 
all its factors, and is therefore given by 6(B) = 60(B) + 61(B) + • • • + 6k(B), where 6i(B) is the number of 
times that Xi appears in B. A Lie bracket B is said to be bad if 60(B) is odd and 6t(B) is even, i £ {1, . . . , k}. 
Otherwise, B is said to be good. To make precise sense of these notions (degree, bad, good) one must resort 
to the concept of free Lie algebras, but it should be clear from the context what we mean here (see [32] for 
a detailed discussion). 

Theorem 4.5. Assume that the system (|4.1|l is locally accessible from m £ M . If every bad Lie bracket B in 
{<7, 771 , . . . , rjk} evaluated at m can be put as an R-linear combination of good Lie brackets in {a, r)i, . . . , rjk} 
of lower degree than B and elements in the kernel of the anchor map p at m, kerp(m), then the system is 
locally controllable from m. 

Proof. Under the hypothesis of the theorem, any bad Lie bracket of the vector fields {/ = p(<r),gi = 
p(f]i), . . . ,gk — p(Vk)} evaluated at m can be written a K-linear combination of good Lie brackets in 
{/, g-y, . . . ,gk} of lower degree. The application of Theorem 7.3 in gives the result. □ 

4.1 Control systems related by a morphism of Lie algebroids 

Let r : E — > M and r : E — > M be two Lie algebroids. Consider a control system defined on E with drift 
section a and independent input sections rji, i = 1, . . . , fc, and a control system defined on E with drift section 
a 7 and independent input sections rjj, j = 1, . . . , fc. Let "J : E — > E be a morphism of Lie algebroids. Let 
spanjr^ , . . . rf^\ denote the linear subbundle of E generated by the sections rjj, j = 1, . . . , fc, and consider the 
affine subbundle £ + spa,n{rj 1 , . . . rf^}. 

Definition 4.6. The control problems on E and E are weakly ^-related (or ^ maps the system on E onto 
the system on E) if ^>(a(m)) 6 £ (ip(m)) + spanj^, . . . rj- k -}(4'( m )) an d ^(Vi( m )) £ s P an {?7ij ■ ■ ■'7fe}('0( TO )) f or 
all i — 1, . . . , fc and all m s M . 

Equivalently, two systems are weakly ^-related if there exist functions Cj and V on M such that 

k k 

* o a = £ o ip + V" V(Tfj o 1/)) and * o r)i = V] Cf '(rfj o tp) , i = l,...,k. (4.3) 

j=i 3=1 

Definition 4.7. The control problems on E and E are ^-related if a is ^-related to a section of E with 
values in £ + spanf?^, . . .rjj} and Tji is ^> -related to a section of E with values in span-fr^, . . ■fh-}, for all 
i = 1, . . . , fc. 

Equivalently, two systems are ^-related if there exist functions C* and V on M such that 

k k 

* o a = £ o -0 + y^ffirjj) o -0 and $ot)j =V(Ci^)of i = l,...,fc. (4.4) 

i=i 3=1 

Clearly, ^-related systems are also weakly ^-related. The following result establishes the relation between 
the controllability properties of weakly ^-related systems. 

Proposition 4.8. Let ^ : E — > E be a morphism of Lie algebroids such that the associated base map tp is open. 
Consider two control systems on E and E that are weakly ^> -related. If the system on E is locally accessible 
(respectively locally controllable) from m £ M, then the system on E is locally accessible (respectively locally 
controllable) from ip(m). 
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Proof. First, note that if two control systems are ^-related, then the image by the morphism ^ of a solution 
a(t) of the system on E with control functions 1 1— > Ui(t), i = 1, . . . , k is a solution of the system on E with 
control functions t h- > Uj(t) = lP(r(a(t))) + J2i=i ( T ( a (t))) u i J = 1, - - • , fe- 

Let m G M and take F an open neighborhood of ^>(m) in M. Consider the open neighborhood V — 
ifi^lV) of m. If the system on E is locally accessible from m, then there is a non-empty open set contained 
in 0^(m, < T) (where Jl^j (m, < T) denotes the reachable set in M starting from m in time less than or 
equal to T, see the appendix). Since if) is an open map and ^(IR]j^(m, t)) C CR^- (^(m), t), it follows that 

= ifi{0) C X^(ip(m), < T), and thus the syste m on E is locally accessible (respectively controllable) from 
ip(m). The argument for the controllable case is analogous. □ 

A interesting particular case occurs when \I/ is an isomorphism between the fibers of the Lie algebroids. 
In such a case, if two systems are ^/-related, one can see that the sufficient conditions for local accessibility 
(respectively controllability) are either simultaneously satisfied on E and E or simultaneously not satisfied. 

Proposition 4.9. Let ^ : E — > E be a morphism of Lie algebroids which is an isomorphism on each fiber. 
Consider two control systems on E and E, with k > k, that are ^-related. Let m £ M . Then 

(a) Lie({a-,77i, . . .,rj k })(m) = E m if and only if Lie({a, rj x , . . -,%})(^(m)) = E^ {m) , 

(b) Every bad Lie bracket in {cr, rji, . . . ,r]k} evaluated at m can be put as an M.-linear combination of good 
Lie brackets of lower degree if and only if every bad Lie bracket in {a, rjj} evaluated at ip{m) £ M 
can be put as an M.-linear combination of good Lie brackets of lower degree. 

Proof. Since the control systems are ^-related and \& is an isomorphism on each fiber, we have that k = k, 
and one can assume without loss of generality that the drift sections a and a and the input sections r\i and 
fj i satisfy 

^ o a = er o tp and \& o rji = fj i o ip , i = 1, . . . , fe. (4-5) 

From the properties of a morphism of Lie algebroids (see Section ITTJl . we deduce that the Lie bracket of any 
subset of sections in {cr, r/i, . . . , r/k} is \l/-related to the Lie bracket of the corresponding subset of sections in 
{a, rji, . . . ,rjk\- Using now the fact that both Lie algebroids have fibers of the same dimension, we conclude 
the result. □ 

As an example application, consider the case of a control system on a manifold M invariant under the 
action of a symmetry Lie group G on M. Assume that the action of G on M is free and proper, so that 
7r: M — > M/G = M is a principal fiber bundle. Then the quotient map "J: TM — » TM/G is a morphism 
of Lie algebroids between E = TM — » M and E = TM/G — > T{M/G). Moreover, ^> is an isomorphism 
in every fiber and its associated base map if) is open. Being the control system on E invariant under the 
action of G, it induces a control system on E. From the above results, we conclude that if the reduced 
system satisfies the sufficient conditions for local controllability (respectively accessibility), then the original 
system also satisfies such conditions. Moreover, since the map ip is open, if the reduced system is not locally 
controllable (respectively accessible), then the original system cannot be locally controllable (respectively 
accessible) . 



5 Mechanical control systems 

In this section, we consider control problems defined on the prolongation TE of E and we make use of 
the special geometry of this Lie algebroid to further investigate the controllability properties of the control 
systems defined on it. Given a Lagrangian function L : E — > R on the Lie algebroid, define the associated 
action functional 3 = / t 1 L dt. Consider the following constrained variational problem: find the extremals of 
among the set of admissible curves with fixed endpoints mo and mi in the base M, i.e. curves a : [io, ii] — ► E, 
m(t) = T(a(t)), satisfying 

p(a(t)) — rh(t) , m(to) = mo and m(ti) = mi . 
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One can see that the infinitesimal variations W : [to»*i] — * ^~E of a curve a : [to>^i] ~~ ¥ E corresponding to 
this constrained variational problem are of the form 

W{t) = ^( m (i)K(i)_ + { — + Cfr'toaUt)) ^ 

where a(t) is a curve on E over m(t) which vanishes at t and t\. A simple calculation and the application 
of the Fundamental Theorem of Calculus show that the equations of motion describing the solutions of the 
constrained variational problem are 

d dL 7 pdL _ i dL 

Euler-Lagrange operator. In alternative terms |7l 1241 155 ] . the infinitesimal variations of a curve a : 
[iojii] - > E can be written as the restriction to the curve of the complete lift of a general time-dependent 
section of E. In such a case, the above equations are precisely the components of the Euler-Lagrange operator 
SL : Adm(E) — > E* , which locally reads 

( d dL dL ■ dL\ 

where {e Q } is the dual basis of {e Q }. The equations of motion just read 

5L = 0. 

Equivalently, an admissible curve a : [to,ti] — > E is an extremal of the action functional 3 if the Euler- 
Lagrange operator SL vanishes at the points of the curve in Adm(_B), 1 1— > (a(t), a{i), a(t)). 

Nonholonomic constraints. In the case of a system subject to (linear) nonholonomic constraints, in 
addition to the above data there is also a subbundle D of the Lie algebroid which prescribes the allowed 
velocities for the system. The equations of motion then are given by the application of the Lagrange- 
d'Alembert principle, 

SL e D° , ae D. 

If there is a projector P : E — > D onto D, denoting Q = I — P, the above equations can be rewritten as 

P*(SL) = 0, Q(o)=0. 
Here P* stands for the dual linear map (P*(6),a) = (6, P(a)}, a e E, 6 e E* . 

Control forces. In the presence of external forces, the equations of motion for both the unconstrained and 
the constrained situations have to be modified. Assume that some input forces {9i, . . . , 9 m } C Sec(E*) act 
on the Lagrangian system on E. Then, the equations of motion for the unconstrained control system read 

m 

5L = Y J Ui6 l , (5.1) 

i=i 

and the equations of motion for the nonholonomically constrained control system are 

m 

P*(SL) = ^ U/ P*(0 ; ), Q(a) = 0. (5.2) 
1=1 
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Mechanical control systems. For the remainder of the paper, we will focus our attention on the class of 
Lagrangian control systems (L, {9i, . . . , 9k}) whose Lagrangian function L : E — > K is of the form 

L(a) = — 9(a, a) — V o r(a) , a e E, 

with S : E Xjvf -E — ► R a bundle metric on £ and V a function on M. We denote by {771, ...,%} the input 
sections of £ determined by the control forces {6%, . . . , 9k} via the metric, i.e., r\i = jig (#*-)■ If ^vs denotes 
the spray associated with the Levi-Civita connection V , the controlled equations (|5.1|l can be written as 

k 

a(t) = /(iys (<*(*)) - (grad g V) v (a(t)) + £ m^r,? (a(t)j) . (5.3) 

«=i 

Note that this system is a control problem on the Lie algebroid TE — > P as defined in Sectional This is the 
reason why we will refer to the control problem with data (S, V, {9i, ■ ■ ■ , 9k}) as a mechanical control system 
on a Lie algebroid. Locally, the equations can be written as 

i— 1 

Alternatively, one can describe the dynamical behavior of the mechanical control system by means of an 
equation on E via the covariant derivative. An admissible curve a : t 1— ► a(t) is a solution of the system 15.311 
if and only if 

k 

Vj t) o(t) + grad s V(m(t)) = £ U2 (t)^(m(t)) . (5.4) 

i=l 

Mechanical control systems with constraints. If the mechanical control system (9, V, {9±, . . . , 9k}) 
is subject to the constraints determined by a subbundle D of E, we can do the following. Consider the 
orthogonal decomposition E = D © Z?- 1 , and the associated orthogonal projectors P : E — > D, Q : E — > D . 
Using the fact that 9(P-, •) = 9(-, one can write the controlled equations (|5.2|) as 

k 

^( V a( t) «W) + P (S™ d S V(m(t))) = «i(*)P(t?i(m(t))) , Q(o) = . 

i=l 

In terms of the constrained connection V a r) — P(V~.rf) + V^(Qrj) (cf. Section f2.2|l . we can rewrite this 
equation as V a (t)a(t) + P(grad g V(m(t))) = J2i = iUi(t)P(rji(m(t))), Q(a) = 0. Since the subbundle D is 
geodesically invariant for the connection V, it follows that any integral curve of the spray associated with 
V starting from ao 6 D is entirely contained in D (cf. Definition ^. 8|) . Since the forcing terms in (|5.5|) coming 
from the potential and the inputs belong to D, the same property holds for the total controlled dynamics. 
As a consequence, the controlled equations can be simply stated as 

fc 

V B(t )o(t)+P(gradgV(m(t)))=^u i (t)P(»j i (m(t))), a e D . (5.5) 

»=i 

Note that one can write the controlled dynamics as a control system on the Lie algebroid TE — > E, 

k 

a(t)^p 1 (r^a(t))-P(gT a d s V)^a(t))+Y / ^(mm) v nt))). (5.6) 

2 = 1 

The coordinate expression of these equations is greatly simplified if we take a basis {e Q } = {e a ,eA} of E 
adapted to the orthogonal decomposition E = D © D ± , i.e., D = span{e a }, T) 1 - = spanje^}. Denoting by 
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(y a ) — {y a ,y A ) the induced coordinates, the constraint equations Q(a) = just read y A = 0. The controlled 
equations (|5.5|l are then 



'dx l 



y A = Q. 



Connection control systems. Given that the structure of the controlled equations is the same both in 
the absence l|5.4[l and in the presence of constraints l|5.5[) . we will in general talk about connection control 
systems on t : E — > M. The dynamics of these systems is governed by an equation of the type 

k 

V a(t) a(t) + ij(m(t))) = J2 Ui{t)r)i{m{t))) . (5.7) 

i=l 

Here V is a connection on E, and {77, rji, . . . ,rjk} are sections of E. We will often refer to n as the potential 
energy term in equations (|5.7|) . Associated with this equation, there is always a control system on the Lie 
algebroid TE — > E given by 

k 

d(t) = ^((Tv-^XaW) +X> i (t)» ? r(a(t))). (5-8) 



5.1 Accessibility and controllability notions 

Here we introduce the notions of accessibility and controllability that are specialized to mechanical control 
systems on Lie algebroids. Let m £ M and consider a neighborhood V of m in M. Define the set of reachable 
points in the base manifold M starting from to as 

3t]£f (m, T) = {to' £ M I 3u £ U defined on [0, T] such that the evolution of (|3~5f) 

for a(0) = m satisfies r(a(t)) £ V, t £ [0, T] and r(a(T)) = to'} . 

Alternatively, one may write $J^(to,T) = t(31 t e ^ V \o m ,T)). Denote 

XXt(m,<T) = |J 3&(m,t). 

Definition 5.1. T"/ie system (|5.8|l is locally base accessible from m (respectively, locally base controllable 
from to,) j/[R]^(m, < T) contains a non-empty open set of M (respectively, 3tj^(m, < T) contains a non-empty 
open set of M to which to belongs) for all neighborhoods V of m and all T > 0. If this holds for any to £ M, 
i/ien i/ie system is called locally base accessible (respectively, locally base controllable ). 

In addition to the notions of base accessibility and base controllability, we shall also consider full-state 
accessibility and controllability starting from points of the form m G E, to £ M (note that full-state is 
meant here with regards to E, not to TM). 

Definition 5.2. The system (|5.8|l is locally accessible from to at zero (respectively, locally controllable from 
m at zero,) if (0 m , < T) contains a non-empty open set of E (respectively, 3t^(0 m , < T) contains a non- 
empty open set of E to which m belongs) for all neighborhoods W of m m £7 cm<i aZZ T > 0. J/ i/iis ZioMs 
/or any m £ M, then the system is called locally accessible at zero (respectively, locally controllable at zero). 

The relevance of the above definitions stems from the fact that, frequently, one needs to control a system 
by starting at rest. Nevertheless it is important to notice that not every equilibrium point at to corresponds 
to the point m . Indeed, there might be other relative equilibrium points, explicitly all those points a £ E 
such that p 1 (r(a)) = 0, i.e., a is in the kernel of the anchor map p and F OL {a m ) = 0, a = 1, . . . ,1. 

Finally, we also introduce the notion of accessibility and controllability with regards to a manifold. 
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Definition 5.3. Let ip : M — ► N be an open mapping. The system (|5.8|) is locally base accessible from 
m with regards to N (respectively, locally base controllable from m with regards to N) i/-0(3?]^(m, < T)) 
contains a non-empty open set of N (respectively, ip(!R^i{ m ) — T)) contains a non-empty open set of N to 
which ip(m) belongs) for all neighborhoods V of m and all T > 0. If this holds for any m G M , then the 
system is called locally base accessible with regards to N ( respectively, locally base controllable with regards 
to N). 

Note that base accessibility and controllability with regards to M with Mm : M — > M corresponds to the 
notions of base accessibility and controllability (cf. Definition EH) - Moreover, if the system is base accessible, 
then it is base accessible with regards to N. The analogous implication for base controllability also holds 
true. 

5.2 The structure of the control Lie algebra 

The aim of this section is to show that the analysis of the structure of the control Lie algebra of affme 
connection control systems carried out in |2(J| can be further extended to control systems defined on a Lie 
algebroid. The enabling technical notion exploited here is that of homogeneity. As we will show later, this 
analysis will allow us to enlarge the class of systems to which the accessibility and controllability tests can 
be applied. 

For the purpose of evaluating the brackets of the accessibility subbundle Lie({F — rj v , rj±, . . . , rfj(}) at initial 
states of the form m , m G M, the discussion on the geometry of TE along the zero-section (cf. Section EJJ 
will be most helpful. Since the Lie brackets in the accessibility subbundle of the mechanical control system 
are linear combinations of the brackets of the elements {F, rf[, . . . ,r)¥ ,t]q}, as an intermediate step we will 
first analyze the structure of the subbundle Lie({F, r}\, ... , r]^,r/ v }). 

Let B be a Lie bracket formed with sections of the family X = {F, rf[, . . . , r)%, r/ v }. For each /, consider 
the following sets 

Br'(X) = {B bracket inX | 5(B) = 1} , 
Br ; (X) = {B bracket inX | B G ?;} . 

The notion of primitive bracket will also be useful. Given a bracket B in X, it is clear that we can 
write B = [Bx^B%], with Bi brackets in X. In turn, we can also write B a = [B a i, B a 2\ for a — 1,2, and 
continue these decompositions until we end up with elements belonging to X. The collection of brackets 
Bx, B 2 , Bn, B 12 , . . . are called the components of B. The components of B which do not admit further 
decompositions are called irreducible. A bracket B is called primitive if all of its components are brackets in 
Br_i(X) UBro(X) U {r}. 

We may now recall the following lemma taken from J7j. Although there it is stated for vector fields, the 
proof can be readily extended to sections of a Lie algebroid since it only relies on two facts: (i) the Jacobi 
identity, and (ii) the fact that X has the property that any bracket in Br;(X), I < 2, is identically zero (which 
follows from Lemma I3.3IT 0. 

Lemma 5.4. Any bracket in Bro(X) U Br_i(X) is a finite sum of primitive brackets. 

As a consequence of this lemma, and the fact that all the brackets in Br;(X), with / > 1 vanish when 
evaluated at the zero-section of E, we conclude that the only brackets we need to consider are the primitive 
brackets in Br_i(X) U Bro(X). We do this next. First, observe the computation of the basic brackets <|3.3ll 
and 1)3. 4[l . In particular, notice that [F, a v ] projects to —a. Second, from Lemma EPT iiL we deduce that any 
bracket in Br_i(X) is the vertical lift of a section of E. From Lemma l3~3T iii). we deduce that the brackets B 
belonging to Bro(X) are projectable sections. We will denote by <tb the section to which it projects. Thus, 
we have B(0 m ) = <tb(to), m G M. The following result completely unveils the structure of these brackets. 

Lemma 5.5. Let B G Bro(X) be a primitive bracket. Then either one of the following is true, 

(i) B = [r, Bx] with Bx G Br_i(X). If ax is the section of E such that B\ = a\ , then gb = —o~\. In 
addition, [a 2 , B) = (ax : a 2 ) v for all a 2 G Sec(-E). 
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(ii) B — [Bi,B 2 ] with Bi,B 2 G Br (X). Then, B(O m ) = [<Jb x > < 7 b 2 ] ( m ) / or all m G M. 

Froo/. Let S e Br (X). Then, either B = [r,^] with B 1 6 Br_ x (X) primitive, or B = [Si,B 2 ] with 
Bi,B 2 G Bi'o(X) both primitive. In the first case, eq. I|3.3|l gives B(0 m ) = [T,Bi](O m ) = — oi(m), where 
B\ = a\. From Q3.4|). we also deduce that [crj , -B] = [(T 2 , o"i]] = (oi : o^) 1 '- In the second case, we have 
that -Bi and B 2 are projectable onto ag 1 and <tb 2 respectively, and therefore [Bi,B 2 ] projects to [ctbi, vb 2 ]- 
Consequently B(0 m ) = [B 1 ,B 2 ](0 m ) = [cr Bl ,crB 2 ](m). □ 

Proposition 5.6. Let m G M. Then, 

Lfo({r, tf, • ■ • , vLv v }) n Ver 0m (te) = s^({ v , m ,..., %})H V , 

Lii({r, %WI) n Hor m (TE) = Lfo(Sy^({77, Vi, ■ ■ ■ ,Vk}))(m) . 



Proof. We prove the inclusion D in the first equality by induction. Let us denote Sym ({77, 771 , ...,rjk}) = 
sj>axi{r),rji,...,ri k } and 

Sym W ({77,771, . . . ,?7 fc }) = span{(cri : a 2 ) \ cr t G Sym ( ' l) ({77, 771, . . . ,77/0}) , h + l 2 = 1} . 

The result is trivially true for I — 1, Sym 1 ({77, 771, . . . ,77fe}) v ' C Lie({r, 77^, . . . , n^ , i] v }) . Assume it is true 

for I and let us prove it for I + 1. Take <7j G Sym l *\{r), 771, ... , i]k}), i = 1, 2, with li + l 2 = I + 1 and observe 
that, using (|3.4|) . 

(a 1 :a 2 ) v = [[^,T},^]. 

By induction hypothesis, a\ G Lie({F, 77^, . . . , 77^, 77^}), and therefore, (<ri : a 2 ) v G Lie({r, 77^, . . . , 77^, rfY). 
From the previous discussion, we know that for the opposite inclusion it is sufficient to look at the primitive 
brackets in Br_i(X). Let B = [Bi,B 2 ] G Br_i(X) primitive, with B\ G Br_i(X), B 2 G Bro(X) primitive. 
From Lemma 15. 51 we have either B 2 = [T, B' 2 \ or B 2 — [B' 2 , B 2 ]. In the first case, we have 

B=[B 1 ,[T,B' 2 }}. 

Using again 13. 41 . we conclude that B G Sym({?7,?7i, . . . ,i]k}) v - As for the second case, using the Jacobi 
identity, we have 

B = [Si, [B 2 ,B%]] = —[B' 2 , [B U B' 2 }] + [B> 2 , [B U B^]] . 

Applying repeatedly the above argument to [Bi, B' 2 ] and \B\, B 2 ] until they are expressed in terms of symmet- 
ric products, we see that B can be expressed as a linear combination of elements in Sym({?7, 771, ... , 77^})^, and 
hence we conclude the result. The second equality is a direct consequence of the first one and Lemma l5.5l □ 

Now, consider the set X' — {T — ■ ■ ■ iVk}- As noted before, the elements in Lie(X') are linear 

combinations of the elements in Lic(X). In fact, for each bracket B' of elements in X', let us define the subset 
S(B') C Br(X) formed by all possible brackets B G Br(X) obtained by replacing each occurrence of T — ?Jq 
in B' by either T or 77Q . Then, one can prove by induction (cf. |17j ) that 

B' = J2 (-l) 5fe+l(B) £, (5.9) 

BeS(B') 

where recall that 5k+i(B) stands for the number of occurrences of 77 v in B. Reciprocally, given an element 
B G Br(X), one can determine the bracket B' of elements in X' such that B G S(B') simply by substituting 
each occurrence of T or t)q in B by T — T)q . We denote this operation by pseudoinv(_B) = B' . 

For each k G N, define the following families of sections in E, 

6^.(77; 77!, ... , m ) = {a G Scc(£) I a v = B" , B" = ]T (-ijfc+i^B , B G Br^-^X) primitive} , 

BeS(pseudoinv(B)) 
nBr_i(X)nBr (X) 

d£Ldm VI, ■ ■ ■ , Vk) = W G Sec(E) | a = a B »,B" = ]T (-1) 5 «+^B , B G Br 2fe (X) primitive} . 

Be5(pseudoinv(B)) 
nBr_i(X)nBr (X) 
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Let e ver (7j;r?i, ...,rj k ) = UkenVvelimVi, ■ ■ ->Vk), e ho r(f?;?7i, ...,%) = ^£^^(77; 771, . . . and denote by 
Cver(j?; ill, ■ ■ ■ ,f]k) and Chor^; Vi, ■ ■ ■ iVk), respectively, the subbundles of the Lie algebroid i? generated by 
the latter families. 

Taking into account the previous discussion, we are now ready to compute Lie({T — r] v , rf[ , . . . , 77^ }) for 
a mechanical control system defined on a Lie algebroid. 

Proposition 5.7. Let m £ M. Then, 

Lie"({T - 77", <,...,%"}) n Ver 0m (T£) - C ver (r?; 7ft, ... , 77 fc )(m)^ , 
Lie"({T - tjW, • • • , *7*» n Hor m (T£) - C hor (f7; 77!, ... , % )(m) . 

Proof. From the definition of the families G ver (ji; 771 , . . . , 77^) and Chor(?7; 771, ... , i/fc), one sees that for each 
B 6 Br(X) primitive, we have computed the R-linear combinations from Br(X) that appear along with B 
in the decomposition of pseudoinv(i?) according to (|5.9|l . Since it is only these primitive brackets which 
appear when making Lie brackets of {r — ?7q , 77^, . . . , 77^}, this will generate Lie({r — tjq , 77^, . . . , 77^}) along 
the zero-section of TE. □ 

Remark 5.8. In the absence of potential terms, i.e., 77 = 0, one has that 

C VC r(0;?7i, ...,r] k ) = Sym({77i, . . .,r) k }) , Ci lor (0;?7i, ...,rj k ) = Lie(Sym({?7i, . . .,%})) . 

It is worth noticing that, in this case, C VO r(0; 771, ... , 77^) C Ch or (0; 771, ... , rj k ). This is not true in general. 



5.3 Accessibility and controllability tests 

In this section we merge the notions introduced in Section f5.1l with the results obtained in Section [5. 21 to give 
tests for accessibility and controllability. 

Proposition 5.9. Let m G M and assume the Lie algebroid E is locally transitive at m. Then the mechanical 
control system H5.8(l is 

• locally base accessible from m if Chor(i]', T)i, ... , rjk)(m) + ker p = E m , 

• locally accessible from m at zero if Ch or (?7; Vi, ■ ■ ■ , ?7fc)( m ) + ker /j = E m and C vel (r); 771, ... , r] k )(m) — 
E m . 

Proof. Consider the accessibility subbundle Lie({r — 7/ v , 77^, . . . , 77^ }) in the Lie algebroid TE. Since E is 
locally transitive at m and Ch or (r); 771, ... , r)k)(m) +ker p m = E m by hypothesis, and Chor(??; ■ ■ ■ , ?7fc)( TO ) Q 
Lii({F - r) v ,tf, . . .,»£})(0 TO ) (cf. Proposition E3, then p 1 (Lie({F - V v ,Tff, . . .,<}))(0 m ) = T m (0 M (M)). 
As a consequence, there exists an open connected submanifold N of M containing m such that 0m(N) is an 
integral manifold of p 1 (Lic({r — rj v , 77^, . . . , 77^})). Let A be the maximal integral manifold of E which contains 
0m(N). Since p 1 (Lic({r — r\ v , 77^, . . . , 77^})) is the accessibility distribution, A is invariant under (|5.8|l and the 
system is locally accessible when restricted to A. Thus the set 3£^j(0 m , < T) is open for all U C A neighborhood 
of m and sufficiently small T. Let V be a neighborhood of m in M and define U — r _1 (V) n A. The result 
now follows from the fact that r is an open mapping and hence the set r(3^^(0 m , < T)) C 3^(771, < T) is 
open in M for T sufficiently small. 

As for the second statement, note that the fact that E is locally transitive at m implies that TE is 
locally transitive at m . In addition, since Chor(?7; t]\, . . . , r) k )(m) + kerp = E m and C vcr (rj; rji, ... , rj k )(m) = 
E m by hypothesis, and CWfa 771, . . . ,r) k )(m), C veT (r]; 771, . . . , r] k )(m) v C Lie({r - rf , rft, . . . , V k })(Qm) by 
Proposition 15.71 we conclude Lie({r — r] v ,r]i , . . . ,rj^})(0 m ) + kerp 1 = % m E. The result now follows from 
Theorem PI □ 

In the absence of potential energy terms, the accessibility tests presented above simply read as follows. 

Proposition 5.10. Let m € M and assume the Lie algebroid E is locally transitive at m. A mechanical 
control system (|5.8|) with no potential terms is 
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• locally base accessible from m if Lie(Sym({?7i, . . . , T}k}))(m) + ker p = E, 



• locally accessible from m at zero if Sym({?7i, . . . , r]k})(m) = E m . 

Remark 5.11. Alternatively, if the sufficient condition for locally base accessibility in Proposition l5 . f 01 is not 
met, i.e., Lie(Sym({?7i, . . . , T)k}))(m) + ker p ^ E m , the corresponding proof also yields the following result. 
Let iV denote the maximal integral manifold of Lie(Sym({?7i, . . . ,i]k}))(m) passing through m. Then, for 
each neighborhood V of m in M and each T sufficiently small, !R]Q-(m, < T) C N , and 0l^j(m, < T) contains 
a non-empty open subset of N . 

The notions of good and bad symmetric products can be stated in a similar way as for Lie brackets. We 
say that a symmetric product P in the sections {77,771, ...,77^} is bad if the number of occurrences of each 
77,; in P is even. Otherwise, P is good. Accordingly, (77, : 77,) is bad and ((77 : 77^) : (77, : 77,)) is good. The 
following theorem gives sufficient conditions for local controllability. 

Proposition 5.12. Let m G M. The mechanical control system (|5.8J) is 

• locally base controllable from m if it is locally base accessible from m and every bad symmetric product 
in {77,77!, . . . ,77^.} evaluated at m can be put as an M.-linear combination of good symmetric products of 
lower degree and elements of ker p, 

• locally controllable from m at zero if it is locally accessible from m at zero and every bad symmetric 
product in {77, 771, . . . , 77^.} evaluated at m can be put as an 'M.-linear combination of good symmetric 
products of lower degree. 

Proof. The proof follows from the following considerations. First, note that every bad Lie bracket in {T — 
T] v r}i, . . . ,77^} gives rise to the vertical lift of a bad symmetric product in {77,7/!, . . . ,77^}. Second, observe 
that every good symmetric product in {77,77!, . . . ,77^} can be alternatively written as a good Lie bracket in 
{r — rj v , r\\, . . . , 77^} evaluated at the zero section 0m (modulo a minus sign). The result is now an application 
of Proposition 14. 51 to the setting of mechanical control systems on Lie algebroids. □ 

The corresponding tests for base accessibility and controllability with regards to a manifold can be proved 
in a similar way. 

Proposition 5.13. Let ip : M — ► JV be an open map. Let m 6 M and assume Tp^{p{E m )) = T^^N. Then 
the mechanical control system (|5.8|) is 

• locally base accessible from m with regards to N if Chor(77;7yi, . . . ,r/k)(m) + p~ 1 (kerip*) = E m , 

• locally base controllable from m with regards to N if the system is locally base accessible from m with 
regards to N and every bad symmetric product in {77, 771, . . . , 77^} evaluated at m can be put as an M-linear 
combination of good symmetric products of lower degree and elements of p~ l {kexip 1e ). 

5.4 Mechanical control systems related by a morphism of Lie algebroids 

In this section we study the relation between the controllability properties of two mechanical control systems 
related by a morphism of Lie algebroids. Consider a mechanical control system on r : E — ► M with T — r\ v 
as drift section (where T is a spray) and inputs {77^, . . . , 77^}, and a mechanical control system on r : E — > M 
with r — rj v as drift section (where T is a spray) and inputs {77^, ... ,77^}. 

Let ^ : E — > E be a morphism of Lie algebroids. Assume the two mechanical control systems are weakly 
T\l/-related. Because of homogeneity, one can deduce that TfoT = Fo \P, so that \? maps the corresponding 
associated connection V onto the associated connection V. Moreover, using the definition of morphism of 
Lie algebroids, one can conclude that 

k k 
* o 77 = 77 o * + ^2 b 3 {Vj *) i * Vi = ^2 C i(Vj , i = 1, . . • , fc , 

j=i 3=1 
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for some functions C\ on M, i.e., the relation by TM/ between the vertical lifts of the potential terms and 
the input sections of the mechanical control systems translates into a relation by '5 of the potential terms 
and the input sections themselves. Following the steps of the proof of Proposition ^. 81 one can also infer the 
relationship between the base accessibility and controllability properties of two ^P-weakly related mechanical 
systems. As above, the control functions of the second system are related to the control functions of the first 
one by means of Uj{t) = V(m(t)) + J2i=i ( m (^)) u i{^) ■ 

Proposition 5.14. Let '5 : E — > E be a morphism of Lie algebroids such that the base map tp is an open 
map. Consider two mechanical control systems which are weakly ^-related. If the system on E is locally 
base accessible (respectively locally base controllable) from m then the system on E is locally base accessible 
(respectively locally base controllable) from i/}(m). 

In the particular case when is an isomorphism between the fibers of the Lie algebroids, if the two 
systems are ^-related, then the sufficient conditions for local (base) accessibility (respectively controllability) 
are either simultaneously satisfied on E and E or simultaneously not satisfied, as stated in the following 
result. 

Proposition 5.15. Let \& : E — ► E be a morphism of Lie algebroids which is an isomorphism on each fiber. 
Consider two mechanical control systems on E and E, with k > k, that are ^-related. Let m G M. Then 

(i) C ver (?7;?ji, . . . ,T] k ){m) = E m if and only if C vcr (rj; rj x , . . . ,%)(i/>(m)) = E^ (m) , 

(ii) C hol {r/;r/i, . . . ,ry fe )(m) + kerp = E m if and only if C\ ior (rj;rj 1 , • • ■ >%)(V'(" 1 )) + ker P = E i>(m)> 

(Hi) Every bad symmetric product in {77, 771 , . . . ,r?fc} evaluated at m can be put as an WL-linear combination 
of good symmetric products of lower degree if and only if every bad symmetric product in {rj, rj 1 , . . . , rjj:} 
evaluated at ip(m) can be put as an W-linear combination of good symmetric products of lower degree. 

Proof. Note that the assumptions of the proposition imply k = k. Therefore, the sections r\i and r\ i can be 
chosen to be ^-related 



6 Applications to simple mechanical control systems and semidi- 
rect products 

In this section we show how the formalism of mechanical control systems on Lie algebroids unifies the 
treatment of several situations which have been previously considered in the literature. We recover known 
accessibility and controllability results for simple mechanical control systems and develop some new ones. 

Simple mechanical control systems 

Let Q be a n-dimensional manifold. A simple mechanical control system is defined by a tuple (QjSjVjSF), 
where Q is the manifold of configurations of the system, S is a Riemannian metric on Q (the kinetic energy 
metric of the system), V € C°°(Q) is the potential function and 3 = {F 1 , . . . , F k } is a set of k linearly inde- 
pendent 1-forms on Q, which physically correspond to forces or torques. The dynamics of simple mechanical 
control systems is classically described by the forced Euler-Lagrange's equations 



\& o r)i = r\ i o ip , l<z<fc. 
Taking into account Q5.1l)[) . one can verify that 



(5.10) 



*°C TO (j);f|l, ...,T] k ) = C VO r(?7;7?l, ■ ■ ■ ,Vk) > 

* o C hoi (r);rn, ...,%) = C ixoI (f);rJ 1 , . . . ,%) o if> . 



(5.11a) 
(5.11b) 



The proof now follows from (|5.11|) and the hypothesis that is an isomorphism on every fiber. 



□ 




(6.1) 
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where L : TQ — » K, i(g, <j) = |S(g, g) - V(g) is the La grangian function of the system. 

There are several ways of intrinsically writing these equations. Here, we present a formulation following 
the Lie algebroid formalism explained above. Consider E — TQ, M = Q and the mappings r = tq : TQ — > E, 
p = Id^Q : TQ — > TQ. Then, it is easy to see that TQ is a Lie algebroid with anchor map Idyg. In this 
setting, the forces in J correspond to sections of the dual bundle E* = T*Q. By means of the musical 
isomorphisms associated with the kinetic energy S, we can consider them as sections ^ = {Yy, . . . , Yfc} of 
E = TQ (i.e. vector fields). Then, equations (|6.1|) are equivalently given by 

k 

ait) = p l (T{a{t)) (gradg V) v (a(t)) + £ Ui(t)Y? (a(t)j) , (6.2) 

where T is the second order equation associated with V s . Here, the mapping p 1 is just the identity in TTM. 

The notion of base accessibility and controllability (resp. accessibility and controllability at zero) in M = 
Q precisely corresponds to the concept of configuration accessibility and controllability (resp. accessibility 
and controllability at zero velocity) in Q as introduced in |2l)| . Also, since p is an isomorphism, we conclude 
that the application of Propositions 15.91 and 15.121 to this case just renders the known tests for accessibility 
and controllability |2U| . 

Simple mechanical control systems with symmetry I 

Assume that a simple mechanical control system (Q, S,V,3f) is invariant under the action of a Lie group G 
on Q, 

$ : GxQ — ► Q 

(g,q) I— » <S>{g,q) = <P g (q) = gq. 

Invariance for the control system means that $*9 = 9, &*V — V and <I>*F l = F\ for 1 < i < k and all 
g G G. The orbit through a point q is Orbc(<z) = {gq \ g E G}. We denote by g the Lie algebra of G. For 
any element £ G g, let £q denote the corresponding infinitesimal generator of the group action on Q. Then, 

T q (Orb G (q)) = {£ Q (q) | £ G g} . 

If the action $ is free and proper, we can endow the quotient space Q/G with a manifold structure such that 
the canonical projection ir : Q — ► Q/G is a surjective submersion. Then, we have that Q(Q/G,G,ir) is a 
principal bundle with bundle space Q, base space Q/G, structure group G and projection w. Note that the 
kernel of Tir consists of the vertical tangent vectors, i.e., the vectors tangent to the orbits of G in Q. We 
denote the bundle of vertical vectors by V, with V q — T q (Orbc{q)), q G Q. 

The action $ induces the lifted action of G on TQ, l:Gx TQ — > TQ, defined by $ g = T$ g . Assuming 
$ is free and proper, we have that $ is also free and proper, and hence p : TQ — » TQ/G, p(v q ) — [v q ], is a 
surjective submersion. Consider then the Lie algebroid defined by T = TQ/G, M = Q/G and 

t(K]) = [?]) P(M) = Tn(v q ) . 

It is not difficult to verify that both r and p are well-defined. The Lagrangian function L being G-invariant, 
it induces a reduced Lagrangian function £ : TQ /G — ► M. The invariance of the forces in 5", or equivalently 
the fact that the sections in ^ are invariant, implies that there exist well-defined sections 23 = {Bi : Q/G — > 
TQ/G}, 1 < i < k, such that po 1^ = Hj o 7r. Finally, the invariance of the potential function and the 
Riemannian metric implies that there exists gradg V such that p o gradg = gradg V o tt. In this setting, 
the notion of base accessibility and controllability precisely corresponds to configuration accessibility and 
controllability in Q/G. However, note that the notions of accessibility and controllability at zero in E are 
stronger than the notions of accessibility and controllability at zero velocity in Q/G, since the former ones 
imply that the reachable sets contain open sets in E — TQ/G, whereas the latter only involve open sets in 
T(Q/G). 

If the mechanical control system is defined on a trivial principal fiber bundle Q = G x Q/G, one may 
consider the canonical projection r : G x Q/G — > G. This mapping is open, and the notions of base 
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accessibility and base controllability with regards to G (cf. Section [OJ precisely correspond to the concepts 
of fiber configuration accessibility and fiber configuration controllability as introduced in [§]. Proposition l5.1 
renders appropriate tests to check these properties. 

A different (and natural) question, however, concerns the search for tests of accessibility or controllability 
in Q which make use of the symmetry properties of the mechanical control system. This is what we analyze 
next. 

Simple mechanical control systems with symmetry II 

Now, we will apply the results of Section 15.41 to simple mechanical control systems with symmetry. Denote 
by Ei -> Ah the Lie algebroid TQ -> Q, and by E 2 -> M 2 the Lie algebroid TQ/G -> Q/G. Note that 
both algebroids have fibers of the same dimension n = dimQ. Let \fr : TQ — > TQ/G be the projection 
mapping associated with the lifted action $, "J" = p. It can be easily verified that \& defines a morphism of 
Lie algebroids. The base mapping ip : Q — > Q/G corresponds precisely to the projection tt. Finally, observe 
that is surjective. 

As a consequence of Proposition 15 . 1 5f ii) . we have that the criterion to test base accessibility in Mi = Q 
(cf. Proposition 

C hor (grad g V\V) = Ei(= TQ) , 

is equivalent to 

C hm (grad g V; B) = E 2 (= TQ/G) . 

Hence, in this way we simplify the computational cost to test the controllability properties of the system, 
since one deals with the reduced representation (and therefore, one works in a space of smaller dimension) . 
These are precisely the results obtained in (5J [20] (although here the analysis is more general since we 
consider nontrivial potential terms). The same simplification occurs for accessibility at zero velocity (cf. 
Proposition 15. 15f i11 , 

C V cr (grad g V;M) =TQ if and only if C vcr (gradg V;-B)=TQ/G. 

As for controllability in Q, Proposition 15 . 1 5T iii) ensures that it is enough to check that the bad symmetric 
products in {gradg V, B\, . . . , Bk} are K-linear combinations of good ones in TQ/G. Finally, Proposition l5.14l 
ensures that if the reduced system is not base accessible (resp. controllable), then the original system is not 
base accessible (resp. controllable). 

Semidirect products 

Let g be a real Lie algebra acting transitively on a manifold M, that is, let g — > X(M) be a surjective Lie 
algebra homomorphism mapping each element £ of g to a vector field £m on M. Define the following Lie 
algebroid structure. The bundle space is E = M x g and the mapping r : E — > M is just the projection onto 
the first factor. The anchor map is given by p(m, £) = ^A/( m )- The Lie bracket is defined via the anchor 
map p and the bracket of constant sections. The latter is defined as the constant section corresponding 
to the bracket on g, that is, if o~\(rn) — (to,£i) an d °~2{ m ) = ( m i£2) a^e two constant sections, then 
[o~i,o-2]{m) = (m, [£i,£2]g)- Note that the case of a linear action on a vector space 2J can be treated as a 
semidirect product by taking an orbit of the action as the base manifold M. 

If we identify TE = TM x Tg = TM x g x g using the left multiplication, an element of TE is of the 
form (a,b,v) — ((m, £), (m, rj), (u m ,£, C))- The condition Tr(v) = p(b) simply implies that v m = rjM{ m )- 
Therefore, we can identify TE with M x g x g x g, and the corresponding maps are 

Ti(m, £,,!],() = (m,£), 2r(m,£, 77, () = (m,rj), p 1 (m ) £,ri,C) = (??m(to),£,C) • 

Let (S, V, {9i, . . . , Ok}) be a mechanical control system on the Lie algebroid E. Assume the bundle metric 
S comes from an inner product on g (and therefore does not depend on the base point), S((m, £1), (to, £2)) = 
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S(£i> £2)- For each £ £ g, define adjj : g — > g by S(ad| 771, 772) = S(»7i, [£j ^Ig)- The spray associated with V s 
then reads 

r vs(m-,0 = (m,£,£,ad|£), 
and the controlled equations of motion are explicitly given by 

k 

a — adjj a = - grad g V(m) + u,-77,-(m) . 

»=i 

If one takes constant sections <7j(m.) = (m, £ = 1, 2, the expression of the Levi-Civita connection is 

V| lff2 (m)= (m,i[6,6] a -^(ad+ 1 6 + ad^)) , 
and the symmetric product is 

((Ti : C7 2 )(m) = (m, ~(ad^ £2 + ad| 2 £1)) ■ 

The above-developed tests can be easily applied to this kind of problems in order to determine whether 
or not the system is base accessible (resp. controllable). 

Systems of the above type appears frequently as mechanical systems defined on homogeneous spaces for 
a given group action. If a group G acting transitively on a manifold M, one can consider the Lie algebroid 
TG x M — > G x M, with anchor map p(v g ,m) — (v g ,0 m ). The bracket of the Lie algebroid is just the 
Lie bracket on the manifold G, where the coordinates of M are considered as parameters. Typically, one 
is thinking of mechanical systems defined on TG that depend on certain parameters which are modeled by 
the coordinates of M. These systems are not invariant under the right action on the group G on itself, but 
they are invariant under the action on G on itself and M at the same time. Therefore, one can consider the 
map 'J: TG X M — > Q X M, ^f(v g ,m) = (vgg^ 1 , gm), which is a morphism of Lie algebroids with base map 
ip: G x M — > M, cf>(g,m) = gm. One can verify that the map is an isomorphism in every fiber and that 
its associated base map ip is open. Therefore, one can apply the results of Section 15.41 to determine whether 
a mechanical system on TG depending on certain parameters modeled by M and which is not G-invariant, 
is locally controllable by analyzing the corresponding system on the Lie algebroid q x M. 

An explicit example of the above type of system is the case of a rigid body G = SO (3) subject to control 
forces with fixed direction (in space), which are the elements in M. Then the above procedure precisely 
consists of studying the control problem in body coordinates. The variables in M evolve dynamically due to 
the non-inertial rotating frame. 

7 Conclusions 

We have investigated the controllability properties of control systems defined on Lie algebroids. We have 
established some general controllability results for nonlinear affine control systems. We have also introduced 
the concept of mechanical control system evolving on a Lie algebroid. After defining appropriate accessibility 
and controllability notions, we have investigated sufficient tests guaranteeing them. We have paid special 
attention to the situation where two control systems are related by means of a morphism of Lie algebroids. 
Finally, we have illustrated the results with the class of simple mechanical control systems and the class 
of systems evolving on semidirect products. Future directions of research will include the investigation of 
controllability tests along relative equilibria of mechanical control systems on Lie algebroids and the treatment 
of models that include gyroscopic forces and dissipation. 
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8 Appendix 

Here we gather some basic definitions concerning control systems defined on manifolds 28 . Let Jl\j(m, T) be 
the reachable set from a point mSMat time T > 0, following trajectories which remain in the neighborhood 
V of m in M for t < T. Denote 

3&(m,<T) = |J D&(m,t). 

t<T 

Definition 8.1. The system (|4.1|l is locally accessible from m G M if 3iYf(rn,< T) contains a non-empty 
open set of M for all neighborhoods V of m and all T > 0. If this holds for any m G M, then the system is 
called locally accessible. 

Definition 8.2. The system B4.1JI is locally controllable from m G M ifJl^^m, < T) contains a non-empty 
open set of M to which m belongs for all neighborhoods V of m and all T > 0. If this holds for any m G M , 
then the system is called locally controllable. 



28 



The accessibility algebra G of the control system (|4.2|) is defined as the smallest subalgebra of X(M) 
containing /, <7i, . . . , gfc. It is not difficult to show that every element of 6 is a linear combination of repeated 
Lie brackets of the form 

[X u [X l ^[...,[X 2 ,X 1 ]...]]], 

where Xi £ {f,gi, ■ ■ ■ ,3fc}, 1 < * < I and I £ N. The accessibility distribution C is defined as the distribution 
on M generated by the accessibility algebra C, 

C(m) = span {^(to) | X vector field in 6} , m € M . 
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